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A successful and commonly used ab initio method for the calculation of crystal field levels and magnetic 
anisotropy of lanthanide complexes consists of spin-adapted state-averaged CASSCF calculations followed by 
state interaction with spin-orbit coupling (SI-SO). Based on two observations valid for Ln(III) complexes, 
namely: (i) CASSCF 4f orbitals are expected to change very little when optimized for different states be¬ 
longing to the 4f electronic configuration, (ii) due to strong spin-orbit coupling the total spin is not a good 
quantum number, we show here via a straightforward analysis and direct calculation that the CASSCF/SI-SO 
method can be simplified to a single configuration-averaged HF calculation and one complete active space Cl 
diagonalization, including spin-orbit coupling, on determinant basis. Besides its conceptual simplicity, this 
approach has the advantage that all spin states of the 4P configuration are automatically included in the SO 
coupling, thereby overcoming one of the computational limitations of the existing CASSCF/SI-SO approach. 
As an example, we consider three isostructural complexes [Ln(acac)3(H20)2], Ln = Dy^+, Ho^'*', Er^+, and 
find that the proposed simplified method yields crystal field levels and magnetic g-tensors that are in very 
good agreement with those obtained with CASSCF/SI-SO. 


I. INTRODUCTION 

The ability of single-molecule magnets (SMMs) to dis¬ 
play long-lived spin-polarized states, which are of inter¬ 
est for the development of molecular magnetic memo¬ 
ries, is fundamentally rooted in the details of the elec¬ 
tronic structure of the ground state and the first few 
excited states of open-shell metal complexes.— In par¬ 
ticular, complexes of trivalent lanthanide ions have re¬ 
cently proven to be promising to achieve SMM behav¬ 
ior at higher temperatures than their transition metal 
analogues.— One characteristic feature of the electronic 
structure of lanthanide complexes is that the lowest en¬ 
ergy electronic states do not differ significantly from 
purely ionic states, leading to wave functions that are 
dominated by spin-orbit atomic J-multiplets weakly split 
by the crystal field potential of the surrounding ligands. 
This leads to unquenched 4f orbital angular momentum 
in the ground state, making the magnetism of the open 
shell 4f electrons particularly sensitive to the surround¬ 
ing low-symmetry electrostatic environment via strong 
spin-orbit coupling, therefore often strongly anisotropic. 
The resulting magnetic anisotropy in turn can lead to a 
high spin-reorientation energy barrier in these complexes, 
which is at the origin of slow magnetic relaxation dynam¬ 
ics and SMM behavior— 

Ab initio calculations have proven very useful to help 
unravel the crucial magneto-structural correlations that 
characterize new or potential lanthanide SMMs— The 
only ab initio method currently used for this purpose is a 
combination of the complete active space self consistent 
field method and the state interaction with spin-orbit 
coupling method^ which usually goes by the acronym 
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CASSCF/RASSI-SO in the literature, after its imple¬ 
mentation in the Molcas software package— We shall 
use the shorter name CASSCF/SI-SO in this paper. 

One of the first successful applications of the CASSCF/ 
SI-SO method to the magnetism of Ln(III) complexes 
was the explanation of a non-magnetic ground state 
Kramers doublet in a triangular Dy(III) cluster— S The 
calculations predicted that the local magnetic anisotropy 
of the Dy centers could be described by Ising-type spins 
whose local axes are tangential to the triangle, result¬ 
ing in a net cancellation of the total magnetic mo¬ 
ment. Subsequent studies used angle-dependent mag¬ 
netic susceptibility measurements on single crystals to 
provide direct evidence that such calculations are indeed 
able to predict the direction of the magnetic easy axes 
and associated g-factors in a number of low-symmetry 
Ln(III) complexes— Besides information on mag¬ 
netic anisotropy, these ab initio calculations also pro¬ 
vide crystal field energy levels which can be compared 
with experimental values, accessible through spectro¬ 
scopic techniques— CASSCF/SI-SO calculations have 
also been used to rationalize the direction and extent 
of the magnetic anisotropy in terms of ligand geome¬ 
try and crystal field models based on electrostatic charge 
distributions to guide the design of new lanthanide- 
SMM candidates and to investigate spin and orbital 
magnetization densities in a series of lanthanide sand¬ 
wich complexes 

CASSCF/SI-SO is now widely used to compute spec¬ 
troscopic and magnetic properties of Ln(III) complexes. 
Here we present a critical assessment of this approach 
and propose an alternative which is both a simplifica¬ 
tion and extension. Section HIl reviews how CASSCF/SI- 
SO is applied to Ln(III) complexes. A discussion of 
the characteristic electronic structure of 4f elements sug¬ 
gests that state-dependent orbital flexibility, provided by 
CASSCF, is of minor importance. This is corroborated 
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by an analysis of the common practice of applying state- 
averaged CASSCF to these systems. Based on these 
findings we propose a simplified method based on one 
set of molecular orbitals, obtained from a configuration- 
averaged Hartree-Fock (CAHF) calculation. This ap¬ 
proach allows a great simplification of the subsequent 
state-interaction problem, which can now be formulated 
as a simultaneous diagonalization of Coulomb repulsion 
and spin-orbit coupling in the basis of Slater determi¬ 
nants of the 4P configuration (CASCI-SO). 

In order to test the CAHF/CASCI-SO approach we 
apply it to three example complexes. Some technical 
details of the method are described in Section IIIIl Sec¬ 
tion m compares the results with those obtained with 
CASSCF/SI-SO and it is shown that a very good agree¬ 
ment is found. 


II. CASSCF/SI-SO TREATMENT OF THE 
ELECTRONIC STRUCTURE OF LN(III) COMPLEXES 

We begin by reviewing the CASSCF/SI-SO method 
as it is applied to mononuclear Ln(III) complexes in the 
molecular magnetism literature. (See, for example. Ref.— 
and references therein.) 

The purpose of the CASSCF step is to obtain wave 
functions that can be thought of as corresponding to the 
atomic Russell-Saunders terms, whose degeneracies are 
weakly split by the presence of the ligand environment. 
For a Ln(III) complex whose formal configuration is 4P, 
this is achieved by choosing the active space to consist of 
n electrons in the seven 4f-like orbitals, giving rise to 
Slater determinants. These are spin-adapted into config¬ 
urational state functions (CSF) of definite spin quantum 
numbers S. For each spin manifold a number of CASSCF 
wave functions is then optimized. Spin-orbit coupling 
(SOC) is introduced in the second step (SI-SO) by diag¬ 
onalizing the SOC operator in the basis of the optimized 
CASSCF wave functions. The resulting eigenvectors are 
then used to calculate expectation values of relevant op¬ 
erators, in particular the magnetic moment. 

Due to the importance of SOC in the rare-earth cou¬ 
pling scheme known from atomic theory, according to 
which the strengths of terms in the Hamiltonian are or¬ 
dered as follows: interelectronic repulsion > spin-orbit 
coupling > crystal field potential, one would preferably 
include as many optimized CASSCF spin states as possi¬ 
ble in the SOC diagonalization. The best possible calcu¬ 
lation in this setting would indeed include all spin states 
of the 4D manifold in the spin-orbit mixing, correspond¬ 
ing to what in atomic theory is known as complete inter¬ 
mediate coupling 1 ^ 

Note that CASSCF in general employs different molec¬ 
ular orbitals for different states. Calculating a matrix el¬ 
ement between any two such states can be computation¬ 
ally expensive in large basis sets, due to the mutual non¬ 
orthogonality of the molecular orbitalsi^ The RASSI 
routine^i^ of Molcas deals efficiently with this problem, 


but it is still a computationally demanding task when 
interaction between a large number of CASSCF states 
is required. For example, Dy(HI), the most studied lan¬ 
thanide in single-molecule magnetism, has a 4f® configu¬ 
ration which corresponds to a total of 735 spin states. A 
complete intermediate coupling calculation would require 
interaction between 735 CASSCF states. This has so far 
not been feasible. In fact, in a recent review concerning 
the application of CASSCF/SI-SO to Ln(HI) complexes 
the authors state that, based on experience, current com¬ 
puter capacities limit the number of states to about 300— 

While CASSCF in principle allows individual opti¬ 
mization of Cl roots with respect to molecular orbital 
rotations, this is in practice only feasible for a few of 
the lowest-energy roots of the Cl matrix. When a large 
number of roots is required, the only viable strategy is to 
resort to state-averaged CASSCF, whereby the molecular 
orbitals are optimized to minimize the average energy of 
the required roots. Thus, within the current application 
of the CASSCF/SI-SO strategy to Ln(HI) complexes, 
for each S, a state-averaged CASSCF calculation is per¬ 
formed often including all states with spin S that are 
possible within the dD manifold, giving equal weight to 
all states in the average. It should be noted that the 
molecular orbitals so obtained are completely indepen¬ 
dent from the Cl problem. This can be seen as follows. 
The CASSCF iterative process involves alternating or¬ 
bital rotation and Cl diagonalization steps. The latter 
determines the roots and their energy. But if only the 
average energy of all roots is required, then diagonaliza¬ 
tion is not necessary, because the sum of all eigenvalues is 
always equal to the sum of all diagonal matrix elements 
(i.e. the trace of the Cl matrix): 

dim H5 

^ F;s,^ = TrHs. (1) 

i=l 

A CASSCF optimization of this type is thus mathemati¬ 
cally equivalent to a minimization of the trace of the Cl 
Hamiltonian matrix with respect to molecular orbital ro¬ 
tations. Since the trace is independent of basis choice, 
there is no need to build and diagonalize the Cl matrix 
at every iteration. Indeed, one obtains exactly identical 
results by first performing an SCF minimization of Tr H 5 
and then, using the orbitals so obtained, a single Cl diag¬ 
onalization to obtain state energies and wave functions. 
This shows that the orbital optimization is completely 
decoupled from the Cl problem. 

The reason that this state-averaging procedure works 
well is probably that molecular orbital relaxation be¬ 
tween states of the 4D space is relatively small. In fact, 
the almost pure atomic nature of the 4f valence shell 
would suggest that MCSCF calculations are not required. 
Characteristic for most Ln(HI) complexes is indeed the 
almost complete absence of covalent mixing of 4f atomic 
orbitals with ligand orbitals. It can therefore be argued 
that the strong electron correlation problem is essentially 
an atomic one. In atoms Coulomb repulsion commutes 
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with both and L'^ electronic angular momenta, so 
that if one could start with CSF’s with good quantum 
numbers S and L arising from a given 4f^ configuration, 
i.e. a representation of the Russell-Saunders terms, one 
would expect such basis to be optimal in order to cap¬ 
ture the dominant features of strong electron correlation 
(static correlation) at the atomic level, or in a few cases 
after solving very small diagonalization problems for the 
determination of accurate RusselRSaunders terms as lin¬ 
ear combinations of few LS'-symmetry adapted CSF’s. It 
is therefore arguable that pure atomic spin-orbit multi- 
plets characterized by a total angular momentum quan¬ 
tum number J, obtained by diagonalization of the SOC 
Hamiltonian on the basis of atomic Russell-Saunders 
CSF’s, would represent the most appropriate guess-states 
for subsequent molecular calculations, possibly for sim¬ 
ple SCF calculations to determine the optimal orbitals 
in the presence of the crystal field potential. However, 
current molecular quantum chemistry codes do not work 
that way and one is typically forced to start from a ba¬ 
sis of CSF’s that are only spin-symmetry adapted, thus 
quite far from being atomic states already taking care of 
on-site correlation. 

If this argument is valid, and we are only interested 
in the properties of the lowest spin-orbit multiplet, we 
should expect that in the CASSCF/SI-SO approach 
what captures the relevant electron correlation effects is 
the attempt to reproduce eigenfunctions as closely as 
the rotational symmetry-breaking effect of the ligands al¬ 
lows, pursued via diagonalization of Coulomb repulsion 
in the basis of the CAS same-spin CSF’s. Such attempt 
would arguably be quite independent of orbital optimiza¬ 
tion which is instead crucial to represent the symmetry¬ 
breaking character of the crystal field electrostatic po¬ 
tential, and consequent splitting of the ground atomic 
multiplet. 

The previous discussion naturally suggests to go one 
step further and assume that the averaged 4f orbitals 
will not depend much on the total spin of the wave func¬ 
tion either, so that just one set of orbitals can be used to 
describe the entire 4P manifold of states. The result of 
this simplification is that we can now formulate the Cl 
on a determinantal basis and that the Cl on the electro¬ 
static Hamiltonian and on the SOC, which are separated 
in CASSCF/SI-SO, can now be combined in just one di¬ 
agonalization step. This approach will be described in 
the next section. 


III. SIMPLIFIED APPROACH: CAHF/CASCI-SO 

The method we propose here consists of two steps. In 
the hrst step a set of optimal molecular orbitals is ob¬ 
tained from a suitable configuration-averaged restricted 
HF-type calculation as detailed in subsection IHI HI In the 
second step the optimized orbitals are used to construct 
all Slater determinants of the open 4P shell, regardless of 
which subset was used in the first step to build the aver¬ 


age energy functional. These determinants form the ba¬ 
sis for a generalized configuration-interaction type matrix 
diagonalization which, besides the usual non-relativistic 
Coulomb repulsion operator, also involves the spin-orbit 
coupling operator, as detailed in subsection IHI Cl 

Conhgurational-average methodologies have been ex¬ 
tensively discussed and applied for many years, since the 
work of McWeeny who proposed them to treat excited 
states associated with any number of open shellsj^i^ 
We also note in passing that, as a cheaper alternative 
to CASSCF, complete active space configuration interac¬ 
tion (CASCI) based on molecular orbitals determined in 
a previous step has also been studied by several workers 
in a variety of contexts.— Visser et ah— applied the 
idea of configuration-averaged orbitals to a lanthanide 
crystal-impurity problem in the context of relativistic 
four-components calculations, which therefore contain 
spin-orbit coupling from the very start. Their method 
has never developed, to our knowledge, into a practical 
non-relativistic ab initio approach dedicated to the cal¬ 
culation of crystal field levels and magnetic properties of 
lanthanide complexes. 


A. Hamiltonian 

The Hamiltonian that we use in the present paper is 
identical to that used for CASSCF/SI-SO calculations in 
MolcaSi^ It is given by the second order Douglas-Kroll- 
Hess (DKH) scalar Hamiltonian, combined with the usual 
non-relativistic Coulomb electron repulsion operator and 
an effective one-electron atomic mean-field (AMFI) ap¬ 
proximation of the DKH no-pair spin-orbit operatorA^ 
A detailed description of this Hamiltonian can be found 
in Ref4S. 

We have implemented configuration-averaged Hartree- 
Fock (CAHF) and spin-orbit-inclusive complete-active 
space conhguration interaction (CASCI-SO) modules in 
a local development code, Ceres {Computational Em¬ 
ulator of Rare Earth Systems), which is based on the 
open-shell version— of the Sysmo software— Since in 
Sysmo the integrals of the DKH operators in the atomic 
basis set are not available, all integrals are computed us¬ 
ing the SEWARD program of Molcas 8.0^ and read 
into Ceres to be used by our CAHF and CASCI-SO 
modules, which are described in the next sections. 

B. Configuration-averaged Hartree-Fock (CAHF) orbitals 

In order to obtain a set of molecular orbitals (MO) 
for the subsequent Cl calculation we minimize the aver¬ 
age energy of the states of the 4f* configuration.— There 
are several ways one might choose to do this. For in¬ 
stance, we could build an average-energy functional over 
all Slater determinants with a constant Ms projection of 
the spin angular momentum. Such approach turns out 
to be rather useful, e.g. for debugging purposes. In fact, 
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when Ms = S^ax, the energy functional is equivalent to 
that minimized during a state-averaged CASSCF calcula¬ 
tion performed with as many roots as there are S = ^max 
states within the df* configuration (see Eq. ©), and thus 
the orbitals obtained via this approach should be equiv¬ 
alent to those obtained via a state-averaged CASSCF 
calculation with Molcas. Moreover, for any other Ms 
value our configuration-averaged approach optimizes or¬ 
bitals by mixing two or more spin-manifolds. Clearly, 
such orbitals are not equivalent to any of the state- 
averaged CASSCF optimized orbitals, and it is thus an 
interesting question whether this can lead to significant 
discrepancies between the two approaches. Finally, we 
can average over all Ms spin manifolds and simply obtain 
a fully-averaged SCF problem within the 4f* configura¬ 
tion, in the spirit of the old McWeeny proposal^— which, 
if our reasoning is correct, should yet again lead to nearly 
atomic 4f-like orbitals with no appreciable discrepancies 
from the SA-CASSCF methodology. 

McWeeny has treated the general case of configura¬ 
tional averaging (see, for example, Refj^, §§6.5-6.6). We 
briefly present here the special case of Ms-configuration- 
averaging, and then generalize the results to those re¬ 
ported by McWeeny. We use McWeeny’s density matrix 
notation. 

Let us consider a system with one closed and one open 
shell, having ni spatial orbitals in the closed shell (or 
ni inactive orbitals), n 2 spatial orbitals in the open shell 
(or 712 active orbitals). Let ua be the number of active 
electrons, with Ua spin-up and np spin-down electrons 
(jiA = na + n/s ), thus with fixed Ms = lf2{na — njs) 
value. Averaging the electrostatic Hamiltonian over the 
Slater determinants that can be formed by oc¬ 
cupying the 712 active orbitals with Ua spin-up electrons 
and np spin-down electrons we obtain the energy func- 


tional: 
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-*^av 

1 
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ni n2 ^ 
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(2b) 


712 2 '^2 

u u,v 

(2c) 


where vi = 2, V 2 = ^ are the (average) occupations of 
the two sub-shells, indices i,j run over the inactive space, 
indices u, v run over the active space, and the Coulomb 
{Juv) and exchange {K^v) integrals for the active space 
are weighted by the Mg-dependent coefficients: 




77277 ^ {nA — 1 ) — 2nanp 

n\ (712 - 1) 

71271a (tIA — 1) — 271271q71^ 

77a {^2 - 1 ) 


( 3 ) 


This expression can be rewritten in terms of density 


and integral matrices on the atomic basis 

Aav = 771 Tr[Ri(h + iGi)] + 112 Tr[R2(h + iGa)], ( 4 ) 

where R^ = is the density matrix of shell i, whose 

LCAO coefficients are contained in the columns of T^, h 
is the matrix of the one-electron Hamiltonian, and 

Gi = iiiG(Ri) -|- ii 2 G(R 2 ) 

G2 = 71iG(Ri) + ll2G'^®(R2) 

are Coulomb-exchange matrices, with 

G(R) = J(R)-iK(R) 

G^^(R) = Af®J(R) - A^^K(R), 

and 

J(R)a/3=^R57(7a|^/5) 

(7) 

K(R)a /3 =^R 57 ( 7 a|/ 5 < 5 ) 

^7 


the usual Coulomb and exchange matrices, respectively. 
Greek letters denote basis set functions. 

We are also interested in optimizing the energy func¬ 
tional originally proposed by McWeeny, arising from av¬ 
eraging over the full set of Slater determinants that can 
be obtained in the chosen active space, regardless of the 
spin-projection quantum number Mg. This can be easily 
recovered by multiplying A^® and A^® in Eqs. © by 
the number of Slater determinants („^/ 2 +Ms) La/T-Ms) 
with constant Ms, summing over all possible Ms, and di¬ 
viding by the total number of Slater determinants. This 
leads to: 



2 (77A-1) 
712 (2712 - 1) 


1 (tia - 1 ) 

V2 (2712 - 1 ) 
( 8 ) 


which determines a modified energy functional, given by 
a modified Eq. where the average electron-electron 
repulsion term within the active space (i.e., last term in 
the last line in Eq. ©) is modified using Eq. ®, thus 
becoming: 


' ^ U,V 

The energy functional Eq. @, or equivalently that 
originally proposed by McWeeny that can be obtained 
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from Eq. (0, is akin to the energy functional arising in 
the restricted open-shell HF (ROHF) theoryi^ Thus it 
can be easily shown that a sufficient condition to mini¬ 
mize Eq. m or Eq. O is to build the density matrices 
for inactive and active spaces from the converged self- 
consistent eigenfunctions of an effective Fock-like Hamil¬ 
tonian: 

Feff = aR2FiR2 -I-6R1F2R1 -I-cRs {vi^i — ^2^2) R3 ( 10 ) 

where = h -|- G^, R^ = 1 — R^, and a, b and c are ar¬ 
bitrary real non-zero parameters that can be adjusted to 
improve convergence of the SCF process. This procedure 
has been implemented via a simple modification of the 
open-shell Sysmo code^, by modifying the mean-field 
repulsion potential within the active space by Eqs. © 
or Eqs. (jH]), which enter Eq. m via the G-matrix G 2 
in Eq. ©■ For the time being the only tools we have 
implemented to achieve convergence of the configuration- 
averaged SCF process are a direct inversion of the iter¬ 
ative subspace (DIIS) algorithm,!- and level shifters for 
the open and virtual shells.!! 

C. Spin-orbit-including complete active space 
configuration interaction (CASCI-SO) 

The CAHF orbitals are now used to set up the full 
configuration interaction calculation in the determinantal 
basis of the open shell. Since there are 7 spatial 4f orbitals 
and n electrons distributed among them, the dimension 
of the Cl secular matrix is (^). The largest dimension 
occurs for n = 7 and is 3432. As this is a relatively small 
number, diagonalization of the Cl matrix is fast. 

The Cl routine that we have implemented is based on 
the cr-algorithm of Olsen et ah^ for the scalar (spin- 
conserving) part of the Hamiltonian (i.e., the scalar one- 
electron plus interelectronic Coulomb repulsion terms). 
We also need to include the Cl matrix elements of the 
SOC operator, which is not spin-conserving. 

The AMFI SOC operator can generally be written as 

= ( 11 ) 

i 

where the summation is over all electrons, and t is a 
function of space only. In second quantization this gives 


U^V 

ffi (.^UV ~ '^^uv)^ua^vl 3 ffi (.^UV + )®I/3®t'a] > (12) 

where t„„ are the AMFI integrals, transformed to the 
active molecular orbital basis. The first two terms are 
spin-conserving, and can be handled by the scalar Cl 
algorithm!^ The third and the fourth term introduce 
spin flips and in order to include these we have supple¬ 
mented the original algorithm with routines that handle 
single excitations from a to /3 spin orbitals and vice versa. 


IV. APPLICATION 

This section presents results of calculations on three 
isostructural complexes: [Ln(acac)3(H20)2], Ln = Dy, 
Ho, Er.— We compare the CASSCF/SI-SO method with 
the CAHF/CASCI-SO method. 

We performed single-point calculations on all three 
complexes, using the crystallographic structures4^ ANO- 
RCC basis sets^ were used on all atoms, contracted to 
[9s8p6d4f3g2h] for Dy, Ho, Er, [3s2pld] for C and O, and 
[2slp] for H. 

CASSCF/SI-SO calculations were done with MoLCAS 
8.0i^ The active space consists of the seven Ln 4f or¬ 
bitals, and is occupied by 9 electrons for Dy, 10 elec¬ 
trons for Ho, and 11 electrons for Er. These occupa¬ 
tions correspond to the trivalent oxidation state and give 
rise to a ground spin-orbit multiplet with total angular 
momentum J = 15/2 for Dy(HI), J = 8 for Ho(HI), 
and J = 15/2 for Er(HI). (Within the single Russell- 
Saunders term approximation, these correspond to ®Hi5 
for Dy(HI), % for Ho(HI), and ^I^ for Er(HI)). Crystal 
field splitting of the J multiplets results in eight low- 
lying Kramers doublets (KD’s) for the Dy and Er com¬ 
pounds. The ground J multiplet of the Ho compound, 
being an even-electron system, splits into seventeen non¬ 
degenerate states. State-averaged (SA) CASSCF opti¬ 
mizations were done on the average energy of all states 
belonging to the highest spin, viz. 21 S = 5/2 states for 
Dy, 35 ^ = 2 states for Ho, and 35 S = 3/2 states for 
Er. The resulting wave functions were spin-orbit cou¬ 
pled with the RASSI module of Molcas. We chose not 
to include states of lower spin in the spin-orbit coupling 
calculation for two reasons: First, as mentioned in Sec¬ 
tion m it is computationally not feasible to include all 
spins of the 4f® manifold of Dy(HI) in a RASSI calcula¬ 
tion (this goes for the 4f^° manifold of Ho(HI) as well). 
In practice, one resorts to an approximation, either by 
discarding states above a certain cutoff energy, or by just 
including the highest-spin states only.— Second, this al¬ 
lows to assess the influence of SOC mixing with states of 
lower spin, by comparison with the CASCI-SO results. 

CAHF/CASCI-SO calculations were done following 
the method described in Section mi Two different types 
of HF averaging were considered to obtain the molecular 
orbitals: in the first, averaging was done over all deter¬ 
minants with maximum spin projection Ms, using the 
A coefficients in Eq. ([31). In the second type averaging 
was done over all determinants, using the A coefficients 
in Eq. (jH) . 

Magnetic g-factors were calculated for each Kramers 
doublet of the Dy and Er complexes.— The Ho complex 
is an even-electron system and as such has no Kramers 
doublets. Nevertheless, it is sometimes possible to find 
two quasi-degenerate states and treat them as a pseudo¬ 
doublet, for which g-factors can be calculated. This has 
been done for the two lowest states of the Ho complex 
(Table IHII) . Note that such pseudo-doublets have only 
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TABLE I. [Ln(acac)3(H20)2]: Comparison of total energies, 
shifted by —13328 Hartree for Dy, —13787 Hartree for Ho, 
and —14256 Hartree for Er. GS = ground state energy. 

SA-CASSCF CAHF CAHF 

Ion on Smax on Ms = Smax on all Ms 


Dy 

Eav 

GS 

-1.021652 

-1.057170 

SI-SO 

-1.021652 -0.777037 

GASGI 

-1.057170 -1.055993 

GASGI-SO 


GS 

-1.079553 

-1.086042 

-1.084697 

Ho 

Eav 

-0.477551 

-0.477551 

-0.327231 




GASGI 


GS 

-0.559519 

-0.559519 

-0.558906 



SI-SO 

GASGI-SO 


GS 

-0.589286 

-0.595589 

-0.594871 

Er 

Eav 

-0.866271 

-0.866271 

-0.789330 




GASGI 


GS 

-0.949870 

-0.949870 

-0.949656 



SI-SO 

GASGI-SO 


GS 

-0.982882 

-0.986258 

-0.985998 

one non-zero 

i principal g-factor.— 


Table H] presents calculated average and ground state 


energies, the latter both with and without inclusion of 
SOC. Note that the SOC-free energies in columns 1 and 
2 are identical. This confirms the equivalence of SA- 
CASSCF and CAHF orbitals for the high-spin subspace, 
predicted by Eq. ©• 

The calculated relative energies and magnetic g-factors 
are summarized in Tables inHiYi It is clear that there is 
minimal difference between the results generated by the 
three methods. The largest changes are observed when 
including all spin states in the spin-orbit coupling, as 
opposed to the high-spin states only. Smaller changes 
are observed when using orbitals averaged over all states, 
as opposed to orbitals averaged over the high-spin states 
only. 


V. CONCLUSION 

We have investigated the application of the CAS- 
SCF/SI-SO ab initio method to the calculation of crys¬ 
tal field splitting and magnetic anisotropy in complexes 
of trivalent lanthanide ions. The two main ingredients 
of this method are: (i) Coupling of Slater determinants 
into Russell-Saunders-like terms by configuration inter¬ 
action in the active space df*; (ii) Coupling of those terms 
into J-like multiplets by spin-orbit state interaction. 
CASSCF performs step (i) but uses state-dependent 
molecular orbitals. This complicates step (ii) because 
the SI-SO program has to calculate matrix elements be¬ 
tween states expressed in mutually non-orthogonal or¬ 


TABLE II. [Dy(acac)3(H20)2]: Calculated relative energies 
and g-factors of the Kramers doublets corresponding to the 
crystal-field split J = 15/2 ground multiplet. 


Doublet 

Energy/cm ^ gi 

92 

93 

SA-GASSGF/SI-SO {S = 5/2) 

1 

0.0 

0.01 

0.01 

19.56 

2 

156.4 

0.26 

0.45 

15.70 

3 

234.8 

2.02 

2.87 

11.28 

4 

289.5 

2.27 

5.87 

7.01 

5 

323.3 

2.12 

4.25 

13.84 

6 

417.9 

0.01 

0.13 

16.28 

7 

477.7 

0.04 

0.08 

18.84 

8 

539.7 

0.01 

0.02 

19.22 

GAHF (Ms = 5/2)/GASGI-SO 

1 

0.0 

0.01 

0.01 

19.44 

2 

154.0 

0.26 

0.45 

15.60 

3 

231.9 

1.92 

2.75 

11.22 

4 

285.8 

2.42 

6.02 

6.99 

5 

319.1 

2.05 

4.13 

13.66 

6 

410.8 

0.00 

0.13 

16.21 

7 

468.9 

0.04 

0.07 

18.78 

8 

529.8 

0.01 

0.02 

19.14 

GAHF (all Ms )/GASGI-SO 

1 

0.0 

0.01 

0.01 

19.43 

2 

155.0 

0.25 

0.43 

15.59 

3 

234.2 

1.87 

2.68 

11.23 

4 

288.5 

2.40 

6.05 

6.97 

5 

321.2 

2.05 

4.22 

13.60 

6 

413.3 

0.01 

0.13 

16.23 

7 

471.8 

0.04 

0.07 

18.79 

8 

533.3 

0.01 

0.02 

19.14 


bitals. Based on the fact that “4f molecular orbitals” 
in Ln(III) complexes are almost pure atomic 4f orbitals 
we have suggested that significant state-dependence of 
the CASSCF molecular orbitals is not expected. This is 
corroborated by the already common practice of state¬ 
averaging CASSCF over a large number of 4f" states. If 
so, a single set of 4f-configuration-averaged orbitals may 
be used to represent all states. As a result, steps (i) and 
(ii) may be combined in a convenient single diagonaliza- 
tion on the Slater determinant basis. 
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